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^ ■ Abstract 

o ■ 



We derive the equation of state for hot nuclear matter using Walecka model 
in a nonperturbative formalism. We include here the vacuum polarisation ef- 

Q . 

2 ■ fects arising from the nucleon and scalar mesons through a realignment of the 

vacuum. A ground state structure "with baryon-antibaryon condensates yields 
the results obtained through the relativistic Hartree approximation (RHA) 
of summing baryonic tadpole diagrams. Generalization of such a state to 
include the quantum effects for the scalar meson fields through the c-meson 
condensates amounts to summing over a class of multiloop diagrams. The 
techniques of thermofield dynamics (TFD) method are used for the finite tem- 
perature and finite density calculations. The in-medium nucleon and sigma 
meson masses are also calculated in a self consistent manner. We examine the 
liquid-gas phase transition at lo"w temperatures (~ 20 MeV), as "well as apply 
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the formalism to high temperatures to examine for a possible chiral symmetry 

restoration phase transition. 
PACS number: 21.65.+f,21.30.+y 
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I. INTRODUCTION 



The understanding of hot and dense matter is an interesting and important problem in 
the context of heavy ion colhsion experiments as well as to study different astrophysical 
objects such as neutron stars. Quantum Hadrodynamics (QHD) is a general framework 
which has been extensively used to study the nuclear matter, both at zero |T| and finite 
temperatures 0, as well as to describe the properties of finite nuclei 0. In the Walecka 
model (QHD-I), the nucleons interact via the scalar (a) and vector {uj) mesons with the 
scalar and vector couplings fitted from the saturation density and binding energy of nuclear 
matter The hot nuclear matter has been studied neglecting the Dirac sea 0, i.e., in the 
so called no-sea approximation. 

To study the nuclear matter at zero temperature including the sea effects in the relativis- 
tic Hartree approximation (RHA), one does a self consistent sum of the tadpole diagrams 
for the baryon propagator [Q. There have also been calculations including corrections to 
the binding energy up to two- loops , which are seen to be rather large as compared to the 
one-loop results. However, it is seen that using phenomenological monopole form factors to 
account for the composite nature of the nucleons, such contribution is reduced substantially 
0] so that it is smaller than the one- loop result. However, without inclusion of such form 
factors the mean-field theory is not stable against a perturbative loop expansion. This might 
be because the couplings involved here are too large (of order of 10) and the theory is not 
asymptotically free. Hence nonperturbative techniques need to be developed to consider 
nuclear many-body problems. 

The approximation scheme adopted here is nonperturbative and, it uses a squeezed 
coherent type of construction for the ground state [|^,|^ which amounts to an explicit vacuum 
realignment. The input here is equal-time quantum algebra for the field operators with a 
variational ansatz for the vacuum structure and does not use any perturbative expansion or 
Feynman diagrams. We have earlier seen that this correctly yields the results of the Gross- 
Neveu model as obtained by summing an infinite series of one-loop diagrams. It was 
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also seen to reproduce the gap equation in an effective QCD Hamiltonian [jTO[ as obtained 



through the solution of the Schwinger- Dyson equations for the effective quark propagator. 



In an earlier work [jTT|, such a nonperturbative method was applied to consider vacuum 
polarization effects in nuclear matter, where it was shown that a realignment of the ground 
state in nuclear matter with baryon-antibaryon condensates is equivalent to the relativistic 
Hartree approximation (RHA). We had then included the quantum corrections arising from 
the scalar meson in a similar way through sigma condensates, amounting to summimg over 
a class of multiloop diagrams. Recently, the formalism has also been generalised to include 
strange baryons and its effects on the composition of neutron star matter as well as gross 
structural properties of neutron stars |12|. In the present paper, we study hot nuclear 
matter including the vacuum polarisation effects arising from the nucleon and scalar meson 
fields. The method of thermofield dynamics (TFD) is used here to study the "ground state" 
(the state with minimum thermodynamic potential) at finite temperature and density. The 
temperature and density dependent baryon and sigma masses are also calculated in a self- 
consistent manner in the present framework. We note that in the Walecka model, TFD has 
been applied to compute perturbatively quantum corrections to the temperature dependent 



Hartree mean field |]T3[ and effect of such corrections on the equation of state. Here we 
shall however use a nonperturbative variational approach to study symmetric nuclear matter 
within Walecka model using TFD for small temperatures (associated with liquid-gas phase 
transition) as well as at high temperatures to discuss possible chiral restoration transition. 
The ansatz functions involved in such an approach shall be determined through minimisation 
of the thermodynamic potential. 

We organize the paper as follows. In section II, we study the quantum correction effects 
from the nucleon and sigma fields in nuclear matter as simulated through baryon-antibaryon 
and scalar meson condensates. The quantum vacuum in nuclear matter is generalized to 
finite temperatures using the thermofield dynamics method. The corresponding condensate 
functions as well as the functions introduced in the definition of the thermal vacuum are 
determined through the minimisation of the thermodynamic potential. This enables us 



to obtain the equation of state (EOS) and various thermodynamic quantities for the hot 
nuclear matter. In section III, we discuss the results obtained in the present nonperturbative 
framework of inclusion of quantum corrections. Unlike the mean field calculations [0, we 
observe that, at temperatures of around 200 MeV, the decrease in the nucleon mass is slower 
when one includes such quantum correction effects, and need not be indicative of a chiral 
symmetry restoration phase transition. We also study the effect of quantum corrections on 
the liquid gas phase transition. The value of critical temperature for such a phase transition 
is seen to be lowered due to quantum effects arising from the scalar mesons. In section IV, 
we give a brief summary of the present work and discuss possible outlook. 



II. QUANTUM VACUUM IN NUCLEAR MATTER 

We shall start with briefly recapitulating the formalism for studying the nuclear matter 
including the quantum correction effects as arising through a realignment of the vacuum with 



baryon-antibaryon and scalar meson condensates |[TT|]. The Lagrangian density is given as 

+ \mluj^uj^ - -^{d^u, - d,u^){d^^u'' - d'-uj^^). (1) 

In the above, if), a, and cj^ are the fields for the nucleon, a, and uj mesons with masses 
M, TTio-, and respectively. The quartic coupling term in a is necessary for the sigma 



condensates through a vacuum realignment, to exist |Tl[]. We retain the quantum nature of 
both the nucleon and the scalar meson fields, where as, the vector uj- meson is treated as 
a classical field, using the mean field approximation for cu-meson, given as {uo^) = 6^oujq. 
The reason is that without any quartic or any other higher order term for the cj-meson, the 
quantum effects generated due to a;-meson through the present variational ansatz turns out 
to be zero. 

The Hamiltonian density can then be written as 

H = Hn + T^a + T^ij, (2) 



with 



'Hn = ip^-ia ■ V + PM)iP + g^aipip, 



(3a) 



(3b) 



(3c) 



We may now write down the field expansion for the nucleon field ip at time t = as given 
by 

1 



[Ur{k)cir{k) + V,{-k)dU-k)] e^-^dk 



(4) 



(27r)3/2 

with cir and c/^ as the baryon annihilation and antibaryon creation operators with spins r 
and s respectively, and U and V are the spinors associated with the particles and antiparticles 
respectively fill]. Similarly, we may expand the field operator of the scalar field a in terms 
of the creation and annihilation operators, at time t = as 

1 r dk 



(T(a;,0) 



a{k) +a\-k)) e 



ikx 



(5) 



In the above, uj{k) = yk"^ + m^. The perturbative vacuum is defined corresponding to this 
basis through a \ vac) = = c/r | vac) = c\j. \ vac). 

To include the vacuum polarisation effects for hot nuclear matter, we shall now consider 
a trial state with baryon-antibaryon and scalar meson condensates and then generalize the 



same to the finite temperatures and densities [0. We thus explicitly take the ansatz for 
the trial state as 



\F) = UMrlvac), 



(6) 



with 



Up = exp / dk f{k) c\^{k) arsCis{—k) — h.c. 



(7) 



Here = 'w|^(cr ■ k)vis and f{k) is a trial function associated with baryon-antibaryon 
condensates. For the scalar meson sector, = UuUi where Ui = exp(i?| — Bi), [i = I, II). 
Explicitly the Bi are given as 



B\ = j dkJ'^Uk)a\k), Bu^ = y dkg{k)a'\k)a'\-k). (8) 



In the above, a'{k) = Uia{k)Uf^ = a{k) — \j ^i^^/o-(fc) corresponds to a shifted field operator 
associated with the coherent state ||Tl],p[] and satisfies the usual quantum algebra. Further, to 
preserve translational invariance fa{k) has to be proportional to 5{k) and we take fa{k) = 
cro(27r)'^/^(5(fc). o"o corresponds to a classical field of the conventional approach |]TT|. Clearly, 
the ansatz state is not annihilated by the operators, c, and a. However, one can define 
operators, d, d^ and b, related through a Bogoliubov transformation to these operators, 
which will annihilate the state \F). 

We next use the method of thermofield dynamics [l^ to construct the ground state for 



nuclear matter at finite temperature. Here the statistical average of an operator is written as 
an expectation value with respect to a 'thermal vacuum' constructed from operators defined 
on an extended Hilbert space. The 'thermal vacuum' is obtained from the zero temperature 
ground state through a thermal Bogoliubov transformation. We thus generalise the state, 
as given by to finite temperature and density as [p!l|,p| 



\F,P) = U,iP)Uj.mF)- (9) 



The temperature-dependent unitary operators Urj{l3) and UpiP) are given as |]14 



U^{f3) = exp J dke^{k, P)h\k)h\-k) - h.c)j . (10) 

and 

UpiP) =exp(^J dk e4k,(3) dl{k) dl{-k) + e+{k,(3) dir{k)dj,{-k) -h.c^. (11) 

The underlined operators are the operators corresponding to the doubling of the Hilbert 
space that arise in thermofield dynamics method. We shall determine the condensate func- 



tions f{k) and g{k), and the functions 9cr{k,P), 9-{k, (3) and 9+{k,(3) of the thermal vac- 
uum through minimisation of the thermodynamic potential. The thermodynamic potential 
is given as 



Q = -p = e-^S - iipB, 



(12) 



where e and S are the energy- and entropy- densities of the thermal vacuum, and ps is the 
baryon density. On evaluation, the energy density of the thermal vacuum is given as. 



with 



7 



dk 



e = {Ti.)i3 — eisi -\- ^uj ''r 



€{k) cos2/(fe) + ^(Mcos2/(fc) + |fc|sin2/(fc) 



(13) 



Co; = guji^Q 7(27r) ^ J dk{cos^ e+ + sin^ 9-) - ^mlujl, 



(cos^ ^+ -sin^ e_), 
(14a) 



(14b) 



and 



11 r dk 



~ 2 (27r)3 7 ^^^^^^^^S" + cosh25f) a;2(/c)(cosh2gr - sinh2gr) 

+ \mlm + ^\alm + 3A7(/3)2 + \mlal + Xa^, 



cosh29a{k,p) 

(14c) 



with 



m 



The entropy density 



1 r dk 



(27r)3y 2uj{k) 



(cosh25r + sinh25r) cosh29a{k, (3). 



(15) 



S = -7(27r)-^ y cifc[sin2^_ ln(sin2e_) + cos^e_ ln(cos^^_) 
+ sin2^+ ln(sin2^+) + cos2^+ ln(cos2^+)] 

+ {'^T^y^ j dk[cosl?9^\u{cosh^9^) - sinh^^^ ln(sinh^^^)] + 5^ 



(16) 



and the baryon density 
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PB = 7(27r)-3 J dk{cos^ 9+ + sin^ 9.). (17) 

In the above, 7 is the spin isospin degeneracy factor and is equal to 4 for nuclear matter. 
Further, is the contribution to the entropy density from tu-meson. Extremising the 
thermodynamic potential, Q with respect to the condensate function f{k) and the functions 
9:f: corresponding to the nucleon sector yields 

tan2/(fc) = — ^^^^4^ (18) 

and 

sin^ = n\ 7' (19) 

^ exp(/?(e*(A;)T/i*)) + l' ^ ^ 

with e*{k) = [k"^ + M*^)^/^ and fx* = fi — g^^ujQ as the effective energy and effective chemical 
potential, where the effective nucleon mass M* = M + QaO'o. 

For the sigma meson sector, on extremising the thermodynamic potential, the functions 
g{k) and 9^ are obtained as 

tanh2,(.) ^ mP) + ^X^o^ 

and 

sinh^ 9^ = ^^^j;^^; ^'(fc, /?) = {k' + M^m"'- (21) 
In the above, the effective scalar meson mass is given as 

M^{pf = ml + 12XI{f3) + UXal (22) 

with 

1 r dk 1 

^^^^ = J YjeTMjm^ 

It is clear from the equation ( PU]) that in the absence of a quartic coupling no such condensates 
are favoured since the condensate function, g{k) vanishes for A = 0. 
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Then the expression for the energy density becomes 

e = eN + e^ + e<^, (24) 

with 

eN = 7(27r)-3 J dk{e + M*^Y/\sm^ 6^ - cos^ 6+), (25a) 
e^j = gu;UJo l{2n)-^ J ^^(sin^ 9_ + cos^ 9+) - ^mlujl, (25b) 

= \mW + W + / dk{k'' + M,(/5)2)V2 cosh2^^. - ?,\I{(3f. (25c) 

After subtracting out the vacuum contributions (6'±=0, /=0 part) for the nucleon sector, 
one obtains, 

Ae^v = e^„.te + Ae, (26) 

where, 

e/ln^te = 7(27r)-^ / dfc(A;^ + M*^) ^^(sin^ ^„ + sin^ ^^^) (27) 

and 

Ae = -7(27r)-3y dk[{k' + M*^f'^ - {k" + M^f''], (28) 

The contribution arising from the Dirac sea effect given by (|28D , is identical to that 
of summing over the baryonic tadpole diagrams of RHA, before renormalisation. This is 
renormalised by adding the counter terms, as for the zero temperature situation, given as 

4 

ect = E Cn^J, (29) 

n=l 

since finite temperature does not introduce any fresh divergences. 

For the cj-sector, the pure vacuum contribution to is subtracted out, which amounts 
to a normal ordering of the number operator ip'^ip. This yields the usual energy density for 
the w-sector, 
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= Qw^oPT - ^^l^l, (30) 

with 

p^^" = 7(27r)-3 J dk{sm^ 9. - sin^ 9+). (31) 

For the a sector, since there are no additional divergences arising from finite temper- 
atures, we adopt the same renormahzation procedure as in Ref. ||TT|. This yields the gap 
equation for field dependent effective sigma mass, M„{j3), in terms of the renormalised pa- 
rameters as 

M^{f3f = ml + UXna^ + 12Ak/^(M. (/?)), (32) 

where 

Simplifying equation ( |25c| ) and subtracting the vacuum contribution, we obtain the energy 
density for the a, 

Ae. = -mWo + SXn4 + ^ (^In (-|) - - j - 3A,/J 

where M„^o and J/o are the expressions as given by eqs. (^2|) and (|33D with ctq = 0. We 
might note here that the gap equation given by (|32D is identical to that obtained through 
resumming the daisy and superdaisy graphs [|l5l and hence includes higher order corrections 
from the scalar meson field. 

In the absence of the quartic interaction (A^j = 0), equation (0) reduces to 

^e- = l^R^^l (35) 

which refers to the RHA. Also, we note that the sign of Xr must be chosen to be positive, 
because otherwise the energy density would become unbounded from below with vacuum 



fiuctuations |T6|-|I^. 
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One then obtains the renormahsed energy density as 



Cren = C/Site + ^^ren + e<^ + Ae<^, (36) 



with, 



Ae.en = -Y^(M*^ln (^) + M\M - M*) - '^-M\M - M*f 

1 S 9^ 
+ —M{M - M*f - —{M - M*f) (37) 
3 12 

as the contribution from the Dirac sea. The thermodynamic potential, f2, given by equation 
(|l^), after subtracting out the vacuum contributions, is now finite and is given in terms of 
the meson fields, ctq and ujq. 

Extremisation of the thermodynamic potential with respect to the meson fields ctq and 
ujq give the self-consistency conditions for ctq (and hence for the effective nucleon mass, 
M* = M + g^ao), as 

and, for the vector meson field, c^o, as 

9co 7 



J dk{sin^e_ -sin'^e+) (38b) 



ml (27r)3 

where sin^ 9^ as the thermal distribution functions for the baryons and antibaryons, given 
through equation (|T^). 



III. RESULTS AND DISCUSSIONS 

We now proceed with the finite temperature calculations for the nuclear matter. The 
values for = glM^/ml and Cl = glM'^/ml, are given as ^2=183.3,167.5,137.9 and 
(7^=114.7,96.45,63.7, for the RHA and with quantum corrections from sigma meson, for 
Xr=1.8 and 5 respectively [jTT|. These values were fitted from the nuclear matter saturation 
properties as po = 0.193/m~'^ and binding energy as —15.75 MeV. For given values of the 
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temperature and the baryon chemical potential, fi, we calculate the different thermodynamic 
quantities, with the meson fields determined self consistently from equation (p8|). 

We plot the pressure as a function of the baryon density, ps for low temperatures to 
examine the liquid gas phase transition, in figures 1 and 2, for the RHA and A/j=1.8 respec- 
tively. At zero temperature, the pressure decreases with density, reaches a minimum, then 
increases and passes through p = at p = po, where the binding energy per nucleon is a 
minimum. The negative pressure indicates a mechanical instability of the uniform nuclear 
matter of saturation density, po- This has however an interesting physical interpretation. 
For the densities below the saturation density, the uniform nuclear matter is unstable and 
will break up into regions of nuclear matter with density, p = po and zero pressure, sur- 
rounded by regions of vacuum with p = and p = ||19[. This pocket in the pressure versus 
density curve, disappears at the liquid gas phase transition point. For RHA, this transition 
appears to occur at around 23 MeV. Inclusion of sigma condensates reduces this critical 
temperature, to 22 MeV for \r=1.8 and to 21 MeV for \r=5. 

For higher temperatures (beyond the liquid gas phase transition), we plot the EOS for 
RHA and for Xji=1.8 in figures 3 and 4 respectively. These plots show a softening of the EOS 
due to quantum effects arising from the scalar meson sector. Such a softening of the EOS 
with inclusion of the quantum corrections was already observed for the zero temperature 



situation |TT|], which had given rise to a lower value for the incompressibility. For a given ps, 
the pressure has the usual trend of increasing with temperature p[. The pressure for = 
is seen to be already nonzero and appreciable at around a temperature of 150 MeV for RHA 
and 200 MeV for Xji=1.8. This has contributions arising from the thermal distributions of 
baryons and antibaryons, as well as from a nonzero value for the sigma field. 

The magnitude of the scalar meson field, a as obtained through the self-consistency 
condition ( |38a| ), is plotted as a function of the baryon density for various temperatures 
in figure 5 for Ar=1.8. It may be noted that for Pb=0, a field becomes nonzero at a 
temperature of around 160 MeV due to thermal effects, and is observed to have attained an 
appreciable value at 200 MeV due to contributions from the thermal distribution functions. 
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The sigma field attaining a nonzero value was also observed for nuclear matter in the mean 
field approximation in Walecka model [0], which had led to a sharp fall in the effective nucleon 
mass between 150 MeV and 200 MeV. The rapid fall of M'^ with increasing temperature, 
was ascribed to resemble the situation, when the system becomes a dilute gas of baryons in 
a sea of baryon-antibaryon pairs. 

The density and temperature dependent nucleon mass is plotted in figures 6 and 7, for 
RHA and Aij=1.8 respectively. The quantum effects are seen to increase the effective nucleon 
mass. It may be noted that the in-medium baryon mass increases with larger value of the 
quartic coupling. Hence, with inclusion of quantum effects, the rapid fall in effective nucleon 
mass as observed in the mean field calculations is not seen in the present calculation. 
The nucleon mass here becomes about half of its vacuum value at T=250 MeV, as compared 
to M*/M ^ 0.2 in Ref. [0 at T=200 MeV. 

The dependence of the entropy density on the temperature and density for \r=1.8 is 
shown in figure 8. This becomes nonzero for zero baryon density at around a temperature 
of 160 MeV, with contributions from the nonzero value for the sigma field. However, the 
increase with temperature, for p_B=0, is rather gradual and need not be associated with 
a phase transition. A similar behaviour of entropy density becoming nonzero at higher 
temperatures, for p_B=0, was also observed in the Walecka model mean field calculations 



as well as in the QMC model ||20|| . 



Finally, the inclusion of quantum corrections through scalar meson condensates enables 
us to determine the in-medium sigma meson mass in a self-consistent manner, which we 
plot in figure 9 for A/j=1.8. An increase in the quartic coupling, Xr increases the value for 
the effective sigma mass [pJ]) . 

IV. SUMMARY 

To summarize, in the present work, we have studied the hot nuclear matter taking into 
account the vacuum polarisation effects within a nonperturbative variational calculation. 
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The approximation here hes in the ansatz for the ground state. A reahgnment of the ground 
state with baryon-antibaryon condensates takes into account the summing over baryonic 
tadpole diagrams. Generahzation of the vacuum with mesons gives rise to summing over 
multi-loop diagrams. We could have generalized the ansatz to include vector meson (uj)- 
condensates in a similar way. However, such an ansatz with the present Lagrangian would 
lead to a trivial solution for the condensate function (thus giving us the mean field result). 



However, with a quartic term |21| for the u meson, the present ansatz for the ground state 
would give rise to a nontrivial solution for the corresponding gap equation. Including such 
a term, it would be interesting to study the quantum effects arising from the vector meson. 

We would also like to mention that the thermal distribution and the condensate functions 
are obtained here through a minimisation of the thermodynamic potential which included 
interactions with the corresponding effective masses for the cr-meson and nucleon determined 



self-consistently through the equations (0) and ( |38a|) at a given temperature. 

We have looked into the liquid-gas phase transition, and observed that the critical tem- 
perature for the phase transition is lowered due to quantum correction effects from the scalar 
meson sector. At high temperatures, with quantum effects, the effective nucleon mass does 
not decrease as rapidly as in the mean field results 0. With increase in temperature, the 
decrease of nucleon mass as well as the increase in the entropy density appear to be rather 
gradual, i.e., without any sudden change and is not indicative of any (chiral or otherwise) 
phase transition. 
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FIGURES 
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FIG. 1. Pressure versus the baryon density for RHA. The disappearance of the pocket at 
temperature of around 23 MeV is indicative of a hquid-gas phase transition. 
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FIG. 2. Pressure versus the baryon density for Xr = 1.8. Quantum corrections due to scalar 
meson sector gives rise to a smaller value of Tc of 22 MeV as compared to that of RHA. 
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FIG. 3. The equation of state for hot nuclear matter for the relativistic Hartree approximation. 
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FIG. 4. The equation of state for = 1.8. The quantum correction from scalar mesons gives 
a softer equation of state. 
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FIG. 5. The meson field strengths as functions of baryon density. 
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FIG. 6. Effective baryon masses in the medium for RHA. 



21 




400 ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — 
0.2 0.4 0.6 0.8 1 

p(fm-3) 



FIG. 7. Effective baryon masses in the medium for Xr = 1.8. The in- medium baryon mass 
increases with quantum correction effects. 
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FIG. 8. Entropy density versus baryon density for A_r = 1.8. 
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